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Abstract 

We show that the OPERA Collaboration have measured the phase 
velocity of the neutrino wave function, based on the analysis of the 
experimental method. On the other hand, the significant advanced 
time displacement 5t reported by OPERA is mainly due to the finite 
life time of the pious (w.r.t. its flying time in the 1 km long tunnel) 
which decay to the muon neutrinos. 

1 . The OPERA Collab. [1] announced to have measured the 'time of 
flight' T {TOFy in [1]) of the neutrino from CERN to Gran Sasso. Dividing 
the distance between these two places L (from BCT to OPERA [1]) by T, 
one can get a 'velocity' whose value turned out to be larger than the speed 
of light in vacuum, c, i.e., L/T > c. Such a result, at first glance, could be a 
challenge to Einstein's Special Relativity, though one is aware that c is one of 
the basic constant in modern physics whose physical meaning is much more 
deeper than light's 'speed' whose physical meaning is quite ambiguous. Only 
in the case of strict classical mechanics, when the position (element of some 
mathematical structure on which differential w.r.t. the parameter 'time' can 
be defined) is considered as the basic observable of a particle, and as a single 
value function of time, the velocity is a well-defined quantity. According to 
Nother's theorem, a basic and well defined space-time related observable is 
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the generator of some kind of space-time symmetry operation. In quantum 
theory, the movement of a particle is described by state vector, and generally 
position is not a good quantum number. So in position representation we en- 
counter a distribution described by wave function of the particle (leaving out 
the more complexity of particle creation and annihilation) . Hence velocity is 
too complex and case-dependent, generally can not be taken as basic observ- 
able related to some kind of space-time symmetry, but as a quantity defined 
by others, depending on the concrete cases. E.g. (suppose Lorentz invari- 
ance kept), for a free particle in energy- momentum eigenstate \E P >, with 
^2 _ p2 _|_ j^2^ its velocity can be defined by P/E, but is not easy to relate 
to some unambiguous distance over time interval as in classical mechanics, 
even with the knowledge of the wave function. Another intuitive example 
is the electron bound in hydrogen atom, whose 3-momentum is not a 'good 
quantum number', i.e., not commuted with the Hamiltonian hence not the 
movement integral. In this case the above definition of velocity fails; and 
because the electron is 'off-shell', if we define some ratio between momentum 
and energy, or space and time, it is very possible to get some 'velocity' which 
is large than c. One may call the electron in such case as tachyon if one likes, 
but we know that the basic physical principles like causality of the hydrogen 
system, e.g., in process of transition between different energy levels, always 
keep. Hence we should pay attention that some 'velocity' defined by ratio 
between some measured space and time intervals is not always straightfor- 
ward pointing to the basic principle about the property of space-time. On 
aware of the above, we can well understand that even the terminology 'speed 
of light in vacuum' is purely a traditional name. Since we should ask what 
kind of 'hght': a photon? a hght pulse? etc. For each case the definition of 
velocity varies. We know that c is further used, e.g., to scale the momentum, 
mass, energy of a free particle so that E"^ = P"^ + M^, etc., etc. However, if 
the wave function of a particle can be deduced from some experiments, one 
of course can discuss some 'displacements' in space as well as time based on 
the wave function. Their ratio could be interesting provided a careful inves- 
tigation on what such kind of measurement the experiment in fact makes, 
and what/whether basic physical principles can be further deduced from it. 

In the following we will recall that the phase speed of the wave function 
can very naturally be larger than c (section 2). With the effort to analyze 
the key experiment method, i.e., fitting the wave forms (average time distri- 
butions in the time interval of the extraction [1]) of the proton as well as the 
neutrino event to get the time displacement w.r.t. 'flying in the speed of c'. 
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we show that the OPERA CoUab. have measured the phase of the particle 
wave function and its speed hence is possibly larger than c (section 3). Since 
neutrino mass is very small, whether phase velocity or group velocity, if prop- 
erly measured, both should not varying from c significantly, no matter larger 
or smaller than c. However, due to the finite life time of the pions (w.r.t. its 
fiying time in the 1 km long tunnel) which decay to the muon neutrinos, the 
wave form for neutrino is deformed from that of the proton. Such an effect 
is estimated to lead to an advanced time displacement of the order of 100ns, 
which coincide with the measured 6t reported in [1] (section 4). 

2. The basic things measured by OPERA CoUab. are two time distri- 
butions in two space-time positions. At position xq (CERN) and time to, 
one measures the proton to deduce the wave function of the neutrino (for the 
problems see section 4) and then at ti and position Xi (Gran Sasso) measures 
the wave function by the OPERA detector (of curse we can only deduce the 
wave function by measuring its mudulus square). To make the discussions 
simple, we do not discuss the neutrino oscillation. The good time correspon- 
dence of the neutrino events and the proton extractions as presented by the 
Collab. [1] indicates that the neutrino is well propagating and could be mod- 
eled as a plane wave for the simplest consideration. We will employ wave 
packet to describe the neutrino in Section 3 and show that the measured 
result L/T is possibly understood as the phase speed of the propagating 
wave function since their key measurement is to calibrate the pace/phase 
for the particle in these two space-time positions. L/T larger than c is the 
straightforward result from that the physical velocity of neutrino is smaller 
than c and well obeys the energy-momentum relation of a massive particle. 
Suppose a likely experiment is done for a much more massive 'dark matter' 
particle, such a phase speed may be significantly larger than c. However, this 
does not to say such measurement is meaningless. Detailed analysis on such 
measurement may help to improve the experiments on neutrino oscillation 
or even absolute value of neutrino mass. 

The wave function is 

and we find that 

^ |Ax| ^ |xi - XqI 
\k\ ~ At ~ ti-to 

This is the well-known phase velocity of the wave of eq. (1). Here one 
encounters the similar problem that Louis de Broglie encountered almost 
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100 years ago [2]. When u and k correspond to energy and momentum of 
a particle, according to the Einstein relation, we will find that this phase 
velocity is never smaller that c but always larger, except for light it equals 
to c. And the wise solution of employing more reasonable group velocity by 
de Broglie is the key corner stone for his matter wave hypothesis which is 
one of the most important turning points from Bohr's old quantum theory 
to the modern quantum mechanics. 

3 . A plane wave as above extends infinitely. If nothing can mark its 
phase, our discussions in Section 2 are just for an 'ideal experiment'. By the 
above we only show that for a plane wave, the interval between two space- 
time positions with the same phase could be space like. For the mechanical 
wave like that on a long shaking rope, we can see the oscillation (say, up 
and down) of the points of the rope with different phases, so that we can 
observe the propagation. It turns out that the time structure of the proton 
PDF (probabihty density function [1], but in fact number density as to be 
clarified) and that of the OPERA neutrino events can mark the phase to be 
measurable, which we will make clear in the following based on the analysis 
of the experiment. It is our key point in this section to make clear what 
the experiment measures with a simple model. Now it is more of reality to 
describe the propagation of neutrino by wave packet. 



where u = u){k) is defined by the on-shell Einstein relation for four momen- 
tum (In this section the space related variables like k, x are not specially 
denoted to show they are S0(3) vectors but indicated). Here we do not dis- 
cuss the details of ^(/c), only this wave packet is assumed well localized: For 
space, even the whole neutrino source system of CERN at Prevessin and the 
detectors in laboratory at Gran Sasso can be taken as point comparing to the 
distance between these two places. So here we need not discuss how well the 
position locahzed which is in fact well investigated by OPERA CoUab. For 
time, we have no problem to say that this particle is localized much smaller 
than 10.5//S, based on the PDF time structure. Same as [1], we employ the 
time when the proton going through the BCT to mark the time of the wave 
packet (e.g., the central peak of the wave packet of Eq. (3)). Taking into 
account the reaction time of the BCT or by carefully analyzing the PDF, 
one can extract a basic time length It (say, around 5 ns according to Fig. 4 
OF [1]). Here a proper ^(/c) need to lead that ^(x,i -|- At) vanishing when 




(3) 
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\At\ > It- Since the wave packet could expand en route of propagation, so 
this requirement is only for the time and place at CERN, where the protons 
employed to produce neutrino by its reaction with the target is carefully 
measured, one of the most important instrument for the this measurement 
is the BCT [1]. 

This is the key point here. We would like to remind the reader that if 
he/she understood that we assign the distribution of the PDF or such kinds 
as the time distribution of the wave (packet) function of each proton as well as 
that of the neutrino, he/she completely misunderstood. Particles produced 
in CERN in the processes described in Section 2 of [1], of which proton is 
the most especially typical, can all be considered as that their wave packet 
size is as small as possible (in time and space) so that the inner structure of 
the wave packet function is not possible seen by all detectors. The PDF time 
structures exactly is the numbers of particles varying with time. However, 
this varying value can mark the PHASE of the wave packet of the particle. 

Since in no way to assume any difference among all the protons, we can 
employ the same form of wave packet to describe all the protons as well as 
neutrinos. Only that, each of the particles (proton/neutrino) corresponding 
to each 'point' {It is the unit) t' at the time axis of PDF distribution has a 
time displacement w.r.t. the beginning of the extraction: 

e-^-^*'*(x,i) ^ J dk *(A;) e-^'^(*+*')+^^^, (4) 

with < t' < lO.b/is. This phase in other word can be called the pace 
of the protons/neutrinos, i.e., to mark their order of 'step' in time of reac- 
tionreaction and reaction of each extraction. (Actually this only marks the 
neutrino's mother particles, mainly pions and kaons, the problems caused by 
the infinite life time during the decay is discussed in details in section 4.) As 
the experiment has put together and makes statistics for all the particles in 
different 'extractions', they are to project on a same t' region by plus/minus 
n time of 50 ms (the time interval between extractions) for their real time 
t. Here we emphasize again that the wave packet is so well located that 
any 'interference' between different wave packets (particles) are neglected, as 
can be assumed by the experiment (or effects as beam interaction has been 
corrected by experiment) . We emphasize this is to point out that those kind 
of effects are of NO relation with our discussions. And we also point out that 
though in experiment each event and wave form has been stamped by the 
real time, but what to be used in analysis is the averaged results obtained by 
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large number of samples of extractions and events projected together on the 
interval of the extraction. We in fact can not say which proton (neutrino) 
has which phase but just employ the different numbers of particles (events) 
in this duration to represent the fact that each wave packet of the particle 
has certain definite phase, marked by t', relative to the, e.g, the beginning 
time of the extraction. 

The wave packet will expand during propagating, but the experiment 
show they still keep some peak structure so that can be detected in Gran 
Sasso at some certain time. This exact value of time but is only used to 
extract the t' to get the event time structure for comparing with the source 
time structure (PDF), then to cahbrate the phase parameter corresponding 
to t'. The method is, as described by the OPERA Collab., to employ the 
maximum likehood method to coinside/tune the same pace for these two 
distributions. This is just to determine which phase group of proton PDF the 
particles/events detected by OPERA belong to. Then it is possible to find the 
same phase positions {to,Xo) and {ti,Xi) when the phase/pace is calibrated 
to coincide. The space-time interval ratio hence can be now understood as 
the velocity of the phase of the wave packet marked by t' moves from (tg, Xq) 
to {ti,xi) (i.e., from CERN to Gran Sasso). xi — xq = L is the length of the 
basehne presented from [1] , while ti — to — T is also possible deduced from 
the measured quantities in [1] . In other words, the experiment is comparing 
the step order of the proton and reaction events, so that to find two space- 
time position with the same phase. Here we would like to emphasize again 
that the PDF and the coinciding neutrino signal distribution is not the form 
of wave function but only show the fact that the wave packet's phase marked 
by t' properly propagates from {xo,to) to {xi,ti). Their interval could be 
space like as discussed in the following. 

The movement of the wave packet from CERN to Gran Sasso, 



in practice is too complex to be employed to calculate the phase velocity 
and information (or 'particle') velocity based on the real form of "^{k), even 
we can get it. Here we only employ the (maybe too) simplified qualitative 
expressions to discuss: Suppose that the effective integral region of k (where 
'^{k) significantly different from zero), Ak, is much smaller than k, i.e.. 



(5) 



Ak « k, which also renders Aoj « u. 



(6) 
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This may be valid in high energies. Now we can rewrite the wave function 
'^{x,t,t') as 

AA;\l/(A;)e~*'^''*^*'-'^*''''^e^*'^'^*^*^'^^. (7) 

Aut' is higher order infinitesimal and neglected from Eq. 7, comparing to a 
long propagation time, since t' is fixed as a small value between to 10.5/is. 
So here t' only appear in the 'carrier' wave. From this quite simple expression 
we see the phase propagates with a velocity ou/k, and the same phase space- 
time position can be marked by t'. While the wave packet propagates at the 
group velocity Au/Ak. The above simplification is just a reproduction of the 
classical discussion from Brillouin' book, chapter 1 [3]. In the most simple 
case the phase velocity is large than c while the group velocity is smaller, 
just to refer Section 2. 

The same narrow k distribution assumption (6) for large k approximately 
leads to 

\<i/{x,t + t')\^^S{vg{t + t')-x), (8) 

with Vg = Au/Ak the group velocity in Eq. (7). This is just the distri- 
bution function of a classical particle and this interpretes why in most high 
energy experiments, one can treat the space-time movement of a single par- 
ticle as classical with velocity P/ E, even with the fact it is almost in energy 
momentum eigenstate rather than position eigenstate. Such a result also con- 
firms what we have emphasized above, that in our treatment, we need not 
to assume anything new from the conventional treatment about the proton 
beam, as well as the pions, kaons, and neutrinos in the production tunnel 
[1] . However, one should be aware that since in quantum mechanics, a single 
microscopic particle is not able to be marked and measured (e.g, the values 
of the position) many times without dramatically changing its state. For a 
static particle fiow, this delta function description is of no difference from 
the plane wave description: no signal of movement like classical one particle 
displacement in space-time can be observed. 

Now comes to the description of the time structure of the proton beam, 
which is also assumed that of the neutrino. This is an inherent property of 
the proton beam system. Here we assumed that many times average (average 
particle number < N{t') >) can cancel the fluctuations from environment, 
i.e., < N{t') > is just the PDF [1]. So the particle system can be described as 
direct product of the the single particle wave functions, with the interactions 
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among them canceled. 

<Nit')> 

which renders 

t + or = 6<''^''>{vg{t + 1') - x). (10) 

Eq. (10) shows that the BCT signal is proportional to the number of the 
protons N{t') as they pass through the BCT at position x and time t+t'. 
When we make the measurement for the times structure of the number of 
particles to mark the same-phase space-time positions as the OPERA exper- 
iment, with the static particle flow rather than one particle movement, we 
are in touch with a space-time distribution B(x,t-|-t'), which is function of 
the phase, i.e., 

<N{t')> 

B = BiYl ^At + t')-kjx) (11) 

The experiments measured two space time positions with the same time 
structure, as one sees the different points of a rope with the same pace of 
oscillation. 

The real information signal, is not the averaged number of < N(t') >, but 
a real time change of this < N[t') >, which propagates with the velocity Vg. 
Supposing a distance large enough so that x/vg » 50/is, then at some time 
we change from the static flow from < Ni{t') >, to < A^2(^') >■ and keep 
the new flow for a time r, with 50 << r << x/vg, then change back to < 
Ni{t') >. These two places can also communicate with light signal, and then 
comparing the propagation time for light and neutrino to get the velocity. 
In fact this is just the analysis based on supernova 1987A observation [4]. 

Such a simple model shows what the measurement the OPERA CoUab. 
have got. Of course, further and deeper analysis, e.g., considering the propa- 
gation as pulse and employ more practical formulations, especially combined 
with oscillation, could help us to improve the measurement of the neutrino 
sector. 

The above discussion in fact is not restricted to the neutrino sector: It is 
a very general discussion based on quantum mechanics principles. However, 
other particle is not possible interact so weak so that can not travel such a 
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long distance. On the other hand, it is not very easy to discuss the causahty 
from the formulae here, since the causal or Feynman propagator is composed 
of positive and negative frequency parts, or say retarded as well as advanced 
parts to keep causality. Here we only employ positive frequency part. So 
a more practical formulation to describe the wave packet of a relativistic 
particle in the quantum field theory framework still need exploration. And a 
precise observable in quantum field theory to define the pace of the particle 
is also in need. 

4- ■ Since the neutrino mass is very small, whether the 'velocity' measured 
by OPERA is the phase one or group one, it should be very near the value of 
c, not as significantly different as reported in [1]. Once the physical meaning 
of the velocity measured by OPERA is clear, the question on what system- 
atics which is missed by the OPERA data analysis cause the advanced time 
displacement 6t ~ 60ns need to be clarified. Such an investigation is very 
important for this kind of measurements is reported to be projected by other 
CoUab., such as MINOS. Here we show that the systematics comes from the 
finite life of the pions which decay to the muon neutrinos, z/^'s. 

As is introduced by [1], the is produced in a 1 km long tunnel, from 
the decay of the production of the proton-target collisions, mainly pion and 
kaon. In hadron-hadron collisions, pion production is dominant. However, 
for higher energy region, the rate of kaon increases. For a neutrino with 
certain energy, the energy of its mother pion or kaon is not quite different in 
average, but the = E/m for kaon is much smaller than that of pion. This 
leads to that the effective (Lorentz retarded) life time of the kaon is much 
smaller than that of the pion and two results are deduced: 1) This effectively 
increase the contribution of kaon to i/^; 2) the kaon can be taken as instantly 
decayed so that the i/^ from kaon can be considered as having the same wave 
form as the proton. 

Let's give the numerical estimation for the 7's and hence the life times of 
various particles which decay to with energy 17 GeV as example, compar- 
ing with their flying time in the 1km tunnel, Tf ~ 0.3 x 10~^. For the case of 
pion, dominantly decaying to fiu^, we choose the case that the momentum of 

in the same direction of that of the pion. We find that 7 must larger than 
several hundreds and its life time T^^ is larger than 10~^s, which is comparable 
with (and larger than) Tf, hence the neutrinos from pions have a deformed 
distribution from that of the proton. For the Kaon sector, with the same ex- 
ample case and same method, we can get the conclusion ~ 10~^s << Tf 
and hence the above conclusions. But the Kg should be treated separately. 
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Since it mainly decays into two pions instantly (comparing to Tf), so its case 
is same as that of the pion. 

Other particles can also contribute, of them muon is very abundant. For 
the case of muon, it decays into z/^ through a 3-body channel. Take the case 
that V/j, energy is half of the mass of muon and momentum is parallel to the 
muon momentum, one can get the 7 is of the order 10^. So 7)^ is of the order 
10~^s >> Tf. So for simplicity and a coarse estimation we can ignore the 
contribution from the muon. If taken into account, the muon number can 
be taken as constant during Tf so will not lead to deformation of the wave 
form. 

So the time of fly in the tunnel Tf is the very key point leads to the 
systematics for the distribution. Mainly the case of pion to z/^ has to be 
considered. Since the neutrino production is proportional to the number 
of the pions, one can propose a exponentially decreasing distribution for the 
neutrino production time in the tunnel. Besides directly to calculate the wave 
form deformation caused by the exponential form of the time distribution of 
the neutrinos from the pion, the effect can be estimated by the displacement 
of the average production time of these neutrinos from the fiat ones in the 
flying interval, with the latter Tf/2. Only considering the effect of pion, it 
leads to the value of the order of 100 ns, which is coincide with the OPERA 
result. The result is calculated from the following equations: 

/o dt t e-*/^- 
dt e-VT. 

Here for simphcity we take Tf as unit. We also emphasize that this effect is 
not sensitive to he wave form of the proton, whether the old long time one 
and the new short time (pulse) one. 

However, here we must point out that such an effect is significant but 
missed by the CoUab. CRUCIALLY depends on the special property of this 
kind of measurement and even more CRUCIALLY depends on the under- 
standing for the physical meaning of this experiment: First, as clarified in 
section 3, this experiment compares phase/pase ('wave form') to measure 
the phase velocity. Second, the recorded neutrino events are quite rare, sta- 
tistically the number of neutrino events recorded by OPERA per proton 
extraction is much much smaller than one. In such a special case, the decay 
process in the tunnel for each extraction is not possible to be tracked by the 
neutrino events (requiring number of neutrino events per extraction > 2), so 
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statistically the phase displacement in the whole tunnel flying will deform 
the recorded neutrino event wave form. Needless to say, if one insists to un- 
derstand this experiment as to directly measure the flying of each neutrino 
recorded by OPERA, in the same way for a classical one, the decay effect 
is very small because 1 — ~ 10~^ and pion only fly not more than 1km, 
which HAS been taken into account in [1]. 

The MINOS [5] experiment is much less affected by the decay effect dis- 
cussed above since its energy is much lower hence a much smaller 7. All the 
mother particles can be taken as to decay instantly. 

5. After the release of the OPERA data, there have been more than 
100 papers discussing this result [6]. This has stimulated the study/review 
of the Lorentz invariance violation in various ways, which could be a good 
window for the future new physics. On the other hand, there are several 
papers devoted to study the systematics of the experiment or suggesting 
the velocity could be unphysical one. However, such kind of investigation 
should include several key points: flrst to clarify what the observable is 
measured by analyzing the basic experimental method; second to see whether 
the measured value agree or not with the suggested, if not, where is the 
systematic. This paper is trying to do in this way. This is a thing that one 
can not escape from for exploring Lorentz invariance violation, not relating to 
any conservative attitude of believing Einstein, but something learned from 
Einstein's attitude to experiments, no matter sooner or later his theories are 
found to be broken through. 



The author thanks Prof. Dr. WANG Meng for explaining some as- 
pects of the OPERA experiments. This work is partially supported by 
NSFC(10935012), SFSD (JQ201101) and SDU (2010JQ006). 

References 

[1] T. Adam et al. [ OPERA Collaboration ], [arXiv:1109.4897 [hep-ex]]. 

[2] Louis de Broglie, Recherches sur la theorie des quanta, PhD Thesis 
(Paris), 1924. 



11 



[3] L. Brillouin, Wave propagation and group velocity, Academic Press, 
1960. 



[4] M. J. Longo, Phys. Rev. D 36 (1987) 3276, Phys. Rev. Lett. 60 (1988) 
173; L. Stodolsky, Phys. Lett. B 201 (1988) 353. The references for 
the observations on the neutrino as well as light signals from Supernova 
1987A are to be found in these three papers. 

[5] P. Adamson et al. [ MiNOS Collaboration ], Phys. Rev. D 76 (2007) 
072005. 

[6] Most of the papers are to be found from 

http://www.slac.stanford.edu/spires/find/hep?c=ARXIV:1109.4897 



12 



This figure "l.jpg" is available in "jpg" format from: 



http://arxiv.0rg/ps/l 1 10.0302v3 



This figure "2.jpg" is available in "jpg" format from: 



http://arxiv.0rg/ps/l 1 10.0302v3 



